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We consider dynamical spin-spin correlation functions in 
the one dimensional repulsive half-filled Hubbard model. We 
propose an exact expression for the two spinon formfactor of 
spin operators. We use this to derive the two spinon contri- 
bution to the dynamical structure factor. 
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I. INTRODUCTION 

The calculation of dynamical correlation functions is 
one of the main challenges in low dimensional theories 
of statistical and condensed matter physics. The large 
distance behaviour (corresponding to very low energies) 
can be effectively extracted for critical theories by using 
bosonization or Bethe Ansatz and conformal field theory 
techniques The results obtained in this way for 

models like the Heisenberg XXZ chain have been success- 
fully applied to neutron scattering experiments of quasi- 
1D magnetic insulators, see e.g. However, from an 

experimental point of view not only very low energies are 
of interest. In magnetic materials described by the spin- 
1/2 Heisenberg model |||| the scattering continuum is 
measured up to energies of twice the exchange, which is 
outside the realm of reliability of the above methods. It 
therefore would be useful to have a method for the cal- 
culation of (Fourier tansforms of) correlation functions 
at all energy scales. Another setting that is outside the 
realm of the above methods are systems with a spec- 
tral gap, see e.g. |^|-^|. Some systems with spectral gap 
can be dealt with by the formfactor approach to quan- 
tum correlation functions 10 13, ij. The advantage of 
this method is that it is in principle not constrained to 
very low energies. The formfactor approach can also be 
implemented for critical systems, in particular this has 
been done for the Heisenberg XXX chain |TJ,|l^]. The 
results are useful for the analysis of the observed two- 
spinon scattering continuum in magnetic insulators like 
KCUF3. In the present paper we consider the problem of 
calculating formfactors in the one dimensional half-filled 
Hubbard model. In a way this model is the "parent" 
of the Heisenberg model used for the description of the 
aforementioned magnetic insulators: the charge sector is 
gapped and the spin sector for large repulsion between 
electrons is described by an effective Heisenberg Hamilto- 
nian. However, if the system is probed at energy scales of 



the order of the Coulomb repulsion e.g. by photoemission 
experiments, which measure the single-particle spectral 
function, the itinerant nature of the system emerges and 
it is necessary to calculate correlation functions in the 
half-filled Hubbard chain in order to describe the exper- 
iments. In the present context it is of particular interest 
to calculate formfactors of the electronic creation and 
annihilation operators, as these results could be directly 
applied to the ARPES data of Jig ], The easier problem 
of determining formfactors of the spin operators in the 
Hubbard model is the subject of the present work. 

We consider the repulsive half filled Hubbard model 
0. The Hamiltonian is (U > 0) 
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This Hamiltonian exhibits an SO(4) symmetry |18|], i.e. 
it commutes with the generators 
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V = ^(-lycj.tCja , rf = - ^(nj, T + n jd - 1) . (3) 
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The complete spectrum of low- lying excitations was de- 
termined in jl9| . It consists of pairs of scattering states of 
four elementary excitations, which form the fundamen- 
tal representation of SU(2) x SU(2). There is a doublet 
of uncharged, gapless spin- 1/2 particles called spinons. 
Their energy and momentum (as functions of the rapid- 
ity variable (3) are given by pti| , |l9t 
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div Jo(w) sin(w2/7/3/7r) 
u> cosh(cj[7) 
du Ji(u>) cos(oj2U (3 / ir) 
uj cosh(uiU) 
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where Jo,i are Bessel functions. The other two elemen- 
tary excitations carry charge ±e but no spin. They are 
called holon and antiholon and have a gap proportional 
to U. Their energy and momentum are 
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e c (k) = 2U + 2 cosfc + 2 



dtv Ji(w)cos(w sinfc)e 
lo cosh(w£7) 



Pc {k) = ir/2 - k - 



duj Jo(w)sin(aj sinfc)e 
a; cosh(a>t7) 



For integrable relativistic quantum field theories one 
generically has a formfactor expansion of the form (^), 
([?]) for correlation functions of local operators. The form- 
(5) factors themselves are determined by the following set of 
axioms ]To| , [Tl[ | 

Axiom 1. The form factors have the symmetry property 



II. DYNAMICAL STRUCTURE FACTOR 

In this paper we consider the dynamical structure fac- 
tor, which is the Fourier transform of the dynamical spin- 
spin correlation function and which is measured by in- 
elastic neutron scattering. The formfactor expansion is 
given by 
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where the labels Sj enumerate the four possible elemen- 
tary excitations. Our notation is the following: e = ± 
denotes a spinon with spin up/down, e = 1,-1 denotes 
an antiholon/holon. As excitations involving holons and 
antiholons have gaps the main contribution to the corre- 
lation function (o) comes from multi-spinon excitations. 
By virtue of related results obtained in other gapless 
models pi]] we expect that the main contribution is due 
to excitations involving only two spinons. Because of the 
spin-5L/(2) symmetry the 2-spinon contribution is of the 
form 
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We propose the following expression for the two-spinon 
formfactor 

/(ft.ft) — = (0|<7+|^ 1 }— = 

c A_(ft - ft)/ [sinh(z7r/4 - ft/2) sinh(^/4 - ft/2)] , 
dt sinh 2 (i[l - 0/iir]) exp(t) 



A-(/3) =exp - 
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sinh(2t) cosh(i) 
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Here c is the usual common constant factor in all form- 
factors. We presently cannot determine its exact numer- 
ical value although it is possible to obtain an estimate by 
considering various sum rules. Let us now provide some 
evidence for the validity of M). 



/(..., ft, Pi+X, ■ ■ ■)e l ,...,e i ,e i + u ...,e n Sel]elXi{Pi ~ ft+l) 

= /(•••, Pi+X, Pi, • ■ ■)e 1 ,...,e' i+1 , E ' i ,...,e n (9) 

Axiom 2. The formfactors fulfil the tensor- valued 
Riemann-Hilbert problem 

/(ft ■ • - ft + 27Ti) ei ... e „ = /(ftft . . ./3„-l) en e I ,..., En _ I 

(10) 

Axiom 3. In the absence of bound states the only singu- 
larities of /(ft, . . . , Pn)si,...,£ n f° r n > 3 are at the points 
Pi — Pj + 7ri, i > j. These singularities are first order 
poles (annihilation poles) with residues 

i Res ./(ft...ft) E1 ,..., £ „ = 

f{Pl, . . . ,Pn-l)e' 1 ,...,e' n _ 2 5e n ,-e„_i • • • <fcj ~ 

sfe'"' 1 - #0 • • ■ C!:^ ~ ft- 2 )J ■ (11) 

Here S^g(P) is the 2-particle scattering matrix of the 
theory under consideration. Although these axioms 
are based on crossing symmetry and relativistic energy- 
momentum relations, there is evidence that the axioms 
hold true even for formfactors of certain nonrelativistc 
lattice models. In particular, it was shown by Pakuliak 
p2| that the known formfactors (T^| of the XXZ Heisen- 
berg lattice model fulfill Axioms 1-3. We believe that this 
will be true for integrable models of statistical mechan- 
ics and condensed matter physics as long as the ground 
state is a singlet under the action of an infinite dimen- 
sional symmetry algebra like a Yangian or U q (sl(2)). In 
particular, the formfactor (||) can be seen to fulfil Axioms 
1-3 



/(/?!,& + 27Ti)__ =f(p 2 , ft)__ , 

/(ft,ft)__ = S__(ft-ft) /(ft, ft). 
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where S (/?) is the dressed S- matrix for scattering of 

spinons in the spin triplet state with S z = — 1 



S—(J3) 



T(-P/2ni)r(l/2 + P/2m) 
T(P/2m) T(l/2- P/2ni) 
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We think that the fact that (g) fulfils Q is a good 
indication for the correctness of the "minimal formfactor" 
A_(P). The full result @ can be checked exactly in the 
limit U — > oo, where the half-filled Hubbard model (|l|) 
reduces to the isotropic spin- 1/2 Heisenberg chain 
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where the exchange is given by J = 1/U. The spinon 
formfactors for this model have been calculated in Jl4],[l5| . 
The dynamical structure factor can be represented in the 
form @, where the sum is over all spinon states. The 
2-spinon contribution is 15 
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x exp(ir7i[p + p(/3i) + p(fh)}) 

x6(lj - e{Pt) - e(/3 2 )) \{0\a+ \foPi}— I 



(15) 



where 



=c' ^_(/3 2 -/?i) 
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Here the dressed energy e(/3) and momentum p((3) are 
given by 



e(j3) = j7r/2cosh/3 , p({3) = arccot(sinh/3) 



(17) 



It is a straightforward to see that in the limit U — » oo 
^(cjjp) indeed reduces to S 2 xxx {'jJ,p)- The £7 — » oo 
limit of e s (/3) and_» ■,(/?) is derived for example in [p5| . 
We believe that (||) is correct for any value of U > 0. 
In p3]] an independent derivation for formfactors of er + 
in the isotropic Heisenberg XXX chain (14) was given 
by using representation theory of a central extension of 
a double of the s/(2)-Yangian. For this purpose <j + can 
be represented as the density of one of the zeroth level 
Yangian generators (Eq in the notation of [^ij ) . Because 
the spin part of the Yangian Q of the Hubbard model 
is the same as in the XXX case, we believe that the XXX 
result can be "lifted" to the half-filled Hubbard model in 
the way presented above. 

In order to perform the integrals over (3\^ in ([?]) we 
need to discuss some properties of the continuum of two- 
spinon excited states. Energy and momentum are given 
by 



£(/3i,A)=e.(/3i) + £«(&) , 
P(/?i,/3 2 ) =p s (/3i)+p s (/3 2 ) mod 2tt, 



(18) 



where (3\^ G (— oo, oo). The upper boundary is obtained 
by taking (3\ — $2 and the lower one by taking (3\ — ±oo. 
We denote the respective dispersion relations by lotj,l(p)- 
In Fig.l we plot the 2-spinon continuum for U = 5. 
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FIG. 1. Continuum of two spinon excitations fur U = 5 . 

In the first Brillouin zone there is a unique solution 
(01,02) to the set of equations 



u> 
P 



PM 



(19) 



as long as u> is chosen within the interval (ul(p), ^u{p))- 
This allows us to taken the integrals over 0x,2 m (@) with 
the result 

S 2 (lo, p ) = c@{w-w L {p)) Q(u>u(p) ~ uj)\A-(0 2 - 0i)\ 2 
x (cosh(/3i/2)cosh(/3 2 /2)) _1 
de s (0 1 ) d Ps (0 2 



dPi d(3 2 
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where Q(x) is the Heaviside function. Due to the compli- 
cated relation between energy/momentum and the spec- 
tral parameter /3 we have not been able to simplify J2(i| ) 
further. It is however in a form that can be readily ana- 
lyzed numerically. Let us first discuss constant momen- 
tum scans, i.e. the behviour of S(to,p) as a function of u> 
for fixed p. At the antiferromagnetic wave vector p = it 
the singularity at zero frequency is 



S(ui, 7r) oc — \/log — for oj — > 

(J V U! 



S(u, 7r) oc y uju(tt) — uj for u> — > wu(jt) 



(21) 



The power-law behaviour in the w->0 limit agrees with 
the result obtained from finite-size corrections and con- 
formal field theory 0. In Fig. 2 we plot S%(u), 7r) on a 

logarithmic scale. One can see that the ^ yj log ^ be- 
haviour holds even for relatively large values of lu. 
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FIG. 2. Constant momentum scan for p = n and U — 5 . 

In Fig. 3 we plot S2 (0.3, p) as a function of the mo- 
mentum p. Obviously the dynamical structure factor is 
nonzero only if intermediate states with energy u> = 0.3 
and momentum p are available i.e. inside the dispersion 
of the 2-spinon continuum shown in Fig.l. 
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discussed in [19] (which have up to a shift in momentum 
by 7r the same dispersion as the two spinon states) . As a 
consequence we can obtain from (||) the two charge- wave 
contribution to pairing and density-density correlation 
functions in the attractive half-filled Hubbard model 

(0|cL, T (*)cLa(*) co 4 (0)co, T (0)|0) , 
(0|n TO (t)no(0)|0) , (23) 

where n m = c* m ^c m .-\ + c* m ^c m .i. The analysis of the 
corresponding Fourier transforms is identical to the one 
carried out above for the dynamical structure factor. 

III. CONCLUSIONS 

We have proposed an exact expression for the two 
spinon formfactor of spin operators in the half- filled Hub- 
bard chain. We explicitly calculated the corresponding 
contribution to the dynamical structure factor. Our pro- 
posal (^J) can be extended to multispinon formfactors of 
spin operators ^7j. We think that very recent devel- 
opments j28| in the calculation of formfactors in inte- 
grable lattice models will make it possible to prove (||) 
rigourously and also lead to the determination of form- 
factors of the fermion operators in the near future. 
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FIG. 3. Constant energy scan for ui = 0.3 and (7 = 5. 

Finally, let us discuss consequences of our results for 
the attractive Hubbard model. It is known that the half- 
filled repulsive Hubbard model can be mapped to the 
half- filled attractive Hubbard model by means of the uni- 
tary transformation defined by 

Uc u tf = c jd , Uc^tf - {-!)* ct T , (22) 

under which H(U) goes to H(—U) and the spin (|J) and 
eta-pairing operators (^|) get interchanged. The repul- 
sive ground state is mapped into the attractive one p6[ , 
which in turn implies that the two spinon states con- 
sidered above get mapped into the charge-wave states 
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